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Abstract. We present a general procedure to renormalize the stress tensor two-point corre- 
lation function on a Minkowski background in position space. The method is shown in detail 
for the case of a free massive scalar field in the standard Minkowski vacuum state, and explicit 
expressions are given for the counterterms and finite parts, which are in full accordance with 
earlier results for the massless case. For the general case in position space, only regularized 
— but not renormalized — results have been obtained previously. After a Fourier transfor- 
mation to momentum space, we also check agreement with a previous calculation there. We 
generalize our results to general Hadamard states. Furthermore, the proposed procedure can 
presumably be generalized to the important case of an inflationary spacetime background, 
where the transition to momentum space is in general not possible. 
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1 Introduction 

The fundamental vertex of the perturbative interaction between any type of matter and 
gravity is given by hahT'^^ , where hah is the graviton field (i.e., the perturbation of a fixed 
background) and Tab is the stress tensor of the matter fields under study. To calculate 
interaction processes between matter and n gravitons it is therefore necessary to obtain n- 
point functions of the stress tensor, and if gravity interacts with a large number of matter 
fields those processes will dominate over graviton self-interactions (this fact can be formalised 
by treating matter fields in an expansion [1, 2]). 

However, because of their extreme smallness quantum effects of gravity are very difficult 
to detect. Except maybe for analogue models of gravity [3], the only field where we can realis- 
tically hope to see any quantum effects is cosmology, where precision measurements became 
recently available [4-6]. Quantum fiuctuations in the early universe are believed to have 
seeded structure formation, and a convenient cosmological observable is the power spectrum 
of primordial gravitational perturbations which manifests itself through temperature fiuctu- 
ations of the cosmological microwave background. Therefore, obtaining quantum corrections 
to this power spectrum is an important task, and a basic building block is the two-point 
function of the stress tensor of matter fields, most importantly of infiaton fields. The general 
case is of course hard, and so in this paper we restrict ourselves to Minkowski spacetime and 
hope to generalize this to curved backgrounds in a later work. On the other hand, the proper 
definition of stress tensor n-point functions is also interesting from a mathematical point of 
view, and having a very concrete procedure may help with a following precise mathematical 
investigation. 

The first obstacle one faces is that the n-point functions obtained in the naive way by 
taking the expectation value of n stress tensor operators are not well defined distributions, 
since they involve products of fields at the same space-time point. One can improve the 
situation by defining a renormalized stress tensor operator which has a finite expectation 
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value, or, alternatively, by working with connected n-point functions (n > 1), which gives a 
finite result as long as all points in the correlation function are not light-like related to each 
other. However, this is still not a well-defined distribution since there remain divergences for 
light-like separations which have to be subtracted by counterterms (p. ex., for the two-point 
function those are terms quadratic in the curvature tensors), and the extraction of those 
divergences and the exhibition of the finite part is in general a hard problem. Only for the 
one-point function, i.e., the expectation value of the stress tensor, the renormalization for a 
generic Hadamard state in a generic curved background is known in full generality [7-11] (this 
can be used to define the above mentioned renormalized stress tensor operator) . For the two- 
point function, only specific examples are known where one can properly separate divergent 
and finite parts and perform renormalization, which include scalar fields on a Minkowski 
background in momentum space [12], massless conformally coupled scalars on conformally 
flat backgrounds [13, 14] and massless minimally coupled scalars in de Sitter space [15], 
while higher n-point functions become tremendously complicated (see [16-20] for examples 
in conformal theories). 

In this paper, we present a procedure to renormalize the two-point correlation function of 
the stress tensor in Minkowski spacetime in a relatively easy way. We explain all steps on the 
concrete example of a massive, minimally coupled scalar field. The method works in position 
space, and so we expect it to be easier generalizable to curved spaces than momentum space 
techniques which are ubiquitous in Minkowski space QFT calculations. In the first section, 
we give a very short summary of Mellin-Barnes integrals which are central to our method, 
and derive some integral formulas which are needed later. ^ The next section gives the first 
step of our procedure, which is the derivation of the (connected) stress tensor two-point 
function for a massive scalar field. This correlation function is regulated in n dimensions and 
given in Mellin-Barnes form. In the second step, done in the section afterwards, differential 
operators are extracted to lower the naive degree of divergence. The resulting expression is 
what in momentum space would be given as the regularized form before renormalization. In 
continuation we split the result into divergent and finite parts (as 71—7-4), and renormalize 
it by adding local counterterms. At the end of the section, we present the fully renormalized 
result. The following section generalizes our results to arbitrary Hadamard states (including 
a concrete simple example), and in the penultimate section the accordance of our results 
with previously obtained special cases is shown. Subsequently, a possible generalization of 
the procedure to curved spacetimes is discussed. 

We employ the "+++" sign convention of [25], and use Latin indices only. The covariant 
derivative with respect to the general metric gab is denoted by Va. In flat space, we also use 
the symbol □ = rj'^^dadi,. We use dimensional regularization in n dimensions. 

2 Integral games 

In this section we explain and derive formulas that will be put to good use in the rest of the 
paper. All of our results are derived by the judicious application of Mellin-Barnes integrals 
— contour integrals in the complex plane which involve products of F functions. An example 
is given by 



'^Mellin transform techniques have also been used for the calculation of Feynman diagrams in momentum 
space (see [21] for an overview), and in the AdS/CFT correspondence (see p. ex. [22-24]). 
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and the corresponding integration path C is shown in figure 1 for a 
c = ^(—2 + i) and d = 0. In general, the integration path goes from 



1 + i) 



4' " 2' 

ioo to ioo, separating 



1 



C 



OV 



Figure 1. The integration path C for the example Mellin-Barnes integral. It separates the complex 
plane in a left (unshaded) and a right (shaded) half such that the left poles (shown in white) all lie 
in the left half and the right poles (shown in black) all lie in the right half. 

left poles (poles of the F functions of the form T[a + z)) from right poles (poles of the T 
functions which are of the form T[a — z)). This is always possible as long as no difference 
between the a's is an integer, so that no left pole coincides with any right pole. Since the 
r functions decay exponentially in imaginary directions, these integrals are well defined, 
and when they decay also in a real direction z — ?■ ±00 they can be evaluated by contour 
integration, summing the resulting series of residues. An important case is given by Barnes' 
Lemma [26, 27] 



V{a + z)T{b + z)T{c- z)T{d- 



dz _ T{a + c)r(a + d)T{h + c)T{h + d) 



IC 27ri T{a + h + c + d) 

Another case which we need in the sequel is given by 

T{a + z)T{-z) ^dz _ 1 
C r(a) ^ 27ri ~ (1 + x)° ' 



(2.2) 



(2.3) 



where the contour can be closed to the left if > 1 and to the right if |x| < 1. 

After this very short introduction to Mellin-Barnes integrals, we use them to calculate 
Fourier transforms that will be important later on. Namely, we need the Fourier transform 
of the Feynman propagator l/(p^ + — ie) in the standard Minkowski vacuum state (with 
the limit e — ?• 0^ understood in the sense of distributions, for which we simply write l/(p^ + 
— iO)). We will start with a general power of the massless case, which is ([21] eq. A. 40, 
[28] eq. 8.715 converted to our conventions) 



d> _.r(f-a) 



1 



(p2 _ iO)« (27r)« 4"7rtF(Q) (x2+i0)i 



(2.4) 
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This is a priori only defined for < < ^, but using the reduction formula 

{x^yp = a{x^y-p (2.5) 

and the fact that the Fourier transform of the d'Alembertian operator is —p^ this can be 
extended to all complex a, except a = and a = f where the T functions diverge and the 
Fourier transform is between a constant and a 6 distribution, 



d"p 

e'P^-^ = 5"(x). (2.6) 
(27r " ^ ^ ^ ^ 



The massive case then follows by application of (2.3) [29] and yields 

e^P^ d> _. f {rr?Y F (f - a - z) r(-z) d^ 



/ 



(p2 + m2 - iO)" (27r)'^ ic (x^ + iO)t-°-^ 4"+^7rtF(a) 27ri 



22-i+°7r2F(a) 



(2.7) 



which for a = 1 is just the Fourier transform of the Feynman propagator [30]. Note that 
in this case the contour may only be closed to the right because in the other half plane the 
integrand is exponentially diverging as z — )• — oo. However, for the following we will only 
need the Mellin-Barnes representation, and use the notation 



2,., 2^.+l-B^(f-l-z)^(-z) dz 



G(x^) = / {m'r{x'Y^'-^ ^2 ' — . (2. 

Jc 4i+^7r2 27ri 

In a similar manner, one obtains the coordinate space Wightman function using 



e(pW + ™')=^'^^^(^^V^-^^^) (2.9) 

27ri + — lU p + ra'^ + lU 



and ([28] eq. 8.730 converted to our conventions)) 

1 n-2a+n+l f + l -i 

^e-^f^ d"a; = \ , , 9(p°)e(V) • (2.10) 

(x2 + i0sgnt)" r(l-f+a)r(a) ^ ^^^(_p2)t-« ^ ^ 

It is just given by G(x2 + i0 sgnt), so that in coordinate space the various two-point functions 
only differ by the type of continuation needed to get around the singularity at x? = 0. 
For derivatives, it follows 

gW(x2) = / (m^)-(x^)-+^-^-t(-l)^ ^ (t - 1 - ^ dz ^2.11) 

For the product of two propagators (and their derivatives), we calculate by shifting the inner 
integration variable 

gW(x2)G«(x2) = I {rr?r{x^r^^-^-'--^^^K{kXz)^^ (2.12) 
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with 

Kik,i,z) = l^r(^'^-i + y-z + k) r(y-z)r(^-i-y + /)r(-y)^ 

_ r{n - 2 + k + 1 - z)r {1 - 1 + k - z) r {'^ - 1 + 1 - z) r{-z) ^^'^^"^ 

~ T{n-2 + k + l-2z) ' 

where we used Barnes' Lemma (2.2) for the evaluation of the Mehin-Barnes integral. It 
is important that in the whole process the contour could be chosen such that no pole was 
traversed. Again, one could evaluate the final integral (2.12) by closing the contour to the 
right (since the integrand diverges exponentially as z — )• — oo), but it is more useful to keep 
the Mellin-Barnes representation. 



3 Regularization 

The action of a free massive and minimally coupled scalar field is well known 

S[g, = / (V>V,,/. + mV') d"x . (3.1) 

By taking a functional derivative with respect to <j), we obtain its equation of motion 

(V^Va - m2) = , (3.2) 
while by deriving with respect to the metric we obtain the stress tensor 

Tab = Va(t>Vb(t> - \gab {V'4>Vc(t> + m^^) , (3-3) 

which, using the equation of motion for (p (3.2), is covariantly conserved V^Tab = 0. This 
can be written as 

Tab{x) = lim Pab{x, vM^My) (3.4) 



Pabix, y) = [ 6l5i - ^-gabg''' V^Vl - -gabn? (3.5) 



with 

and the subscript indicates on which field the derivative acts. 

We now specialize to flat space and consider the connected two-point correlation function 
of the stress tensor 

(T^b{x)T,d{x)) = {0\Tab{x)T,d{x'm - {Tab{x)){TM)) 

= (0| {Tab{x) - (Tabix))) {TM) - {Tcd{x))) |0) ^ - ' 

in some state |0) for the scalar field, which we will take to be the standard Minkowski 
vacuum. The advantage of considering the connected correlation functions is that they are 
finite as long as the points x and x' are not null (light-like) separated, since the divergences 
that arise from 0^(x) (and its derivatives) in the definition of the stress tensor are absent 
in the difference Tab{x) — {Tab{x)). However, in the null separation limit (x — x')^ — )■ the 
correlation functions still diverge after smearing with test functions; they are not well-defined 
distributions in four dimensions. 
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To obtain explicit expressions, we insert the definition of tlie stress tensor (3.3) into tlie 
connected two-point function (3.6) and evaluate the expectation values of k matter fields (f) 
by Wick's formula. To simplify the resulting expression, we take advantage of the fact that 
the operator Pmn{x-, y) (3.5) is symmetric under the interchange of x and y, and it results 

{Tab{x)TM)) = 2 lim Pab{x, y)Pcd{x\ y')G[x, x')G{y, y') , (3.7) 

where 

G{x,y) = {mx)4>{ym- (3.8) 

This two-point function only depends on the coordinate difference x — y because of Poincare 
invariance of the standard Minkowski vacuum |0), and for simplicity we may set y = 0. We 
then get for the two-point function 

{Tab{x)ndm = ^VabVcd [l6{x^G"f + Wx^G'G" + 4(n - 4 + 2m^x^){G'f + m^G^] 

+ 8ria(cVd)b{G'f + ^'^X(a'nb){c)Xd)G'G" + 2,2xaXbXcXd{G"f (^•^) 
- 4 {riabXcXd + XaXbVcd) [4x^{G" f + 4G'G" + m^G'f] , 

where G" = d'^G{x, 0)/9(x^)^ = G"{x'^) and analogue for G' and G. In the above condensed 
notation, the equation of motion for the scalar field (3.2) implies 

ix'^G" + 2nG' = m^G , (3.10) 

and using it the two-point functions may be simplified somewhat; also covariant conservation 
can be checked easily. 

The attentive reader may have noticed that so far the entire discussion has been in terms 
of Wightman functions, while for in-out calculations Feynman propagators (time-ordered 
functions) are needed.^ While they are quite different in momentum space, in coordinate 
space these various functions differ only by the type of analytic continuation needed to 
get around the singularity in the propagator as x^ — )• 0, as shown in the previous section. 
When acting with derivatives on the Feynman (or Dyson) propagator, naively additional 
local terms ~ (5"(x) are obtained, so in this case one should find additional local terms in 
the above expressions for the stress tensor correlation function. However, since a priori the 
time-ordered stress tensor correlation function is not well defined, there is no reason to keep 
those local terms. We will instead assume the working hypothesis that the time-ordered 
two-point function of the stress tensor is obtained by using the Feynman prescription after 
performing the manipulations in this section and the first half of the next section. This is in 
accordance with momentum space calculations, where a derivative da is just multiplication 
with \pa, and the manipulations in the next section just correspond to a reordering of those 
Pa, without additional terms. 

We now put the Mellin-Barnes representation (2.12) for the product of two propagators 
(and their derivatives) into the stress tensor two-point function (3.9). After shifting some 
integration variables, we obtain 

( \rr fnw f( 2x^/ 2^^-n 1 r(n-z)r^(f -z)r(-z) dz 
{Tab{x)T^dm = j^{m ) {X ) r(n + 2 - 2z) ^^'^'^'^2^. ^^'^^^ 

^In the m-m formalism one even needs the full arsenal, including in addition Dyson (anti-time-ordered) 
functions and positive and negative frequency Wightman functions. 
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with 

Tabcdiz) = 2 [ilabVcdin^ - n - 4 - 2z(2n - 1) + Az^) + 'ir]a(cVd)b] {n + l- 2z)(n - 2z) 

- 16^(^^%^(n -z){n- 2z){n + 1 - 2z) 

. VabXcXd + XaXblJcd / n/ n \i- r, n \r ii o\ 

— 4 ^ (n — — 2z)(n — 2 — 2zj(n + 1 — 2zj 

+ g WcX, ^^ + 1 _ _ _ 2^)2 . 

[X ) 

(3.12) 

Although this result looks complicated, it is quite simple and has a very nice property. 
Namely, on physical grounds one expects any stress tensor correlation function to be renor- 
malizable if it fulfils the properties expected from such a function, which is symmetry in 
the index pairs (ab) and (cd) and conservation. No matter what the exact dependence on 
the spacetime coordinates looks like, as long as those properties are fulfilled it is conceivable 
that there is some theory and some quantum state that has this function as the correlation 
function of its stress tensor. It is now easily checked that conservation is independent of z, 
and so it should be possible to rearrange and renormalize the integrand which is simpler than 
the whole (integrated) result. That this in fact can be quite simply done is shown in the next 
section. 

4 Renormalization 

The expression for the two-point function (3.11) derived in the previous section is not yet 
suited for renormalization. The residue of the first pole in the Mellin-Barnes representa- 
tion (3.11) is proportional to (j;^)"", which is more singular than the leading singularity in 
the product of two propagators [G{x'^)]'^ ~ (x^)^""' (2.12). This is of course due to the fact 
that we took derivatives, which in momentum space just correspond to multiplications by pa- 
To reduce the strength of the singularity, we therefore have to extract differential operators 
(which in momentum space corresponds to a reordering of the pa), and to reduce the singu- 
larity to (x^)^"" we see that we need fourth-order differential operators. We are led to the 
exact form of the differential operators we need to extract by the conservation of the stress 
tensor. Of the five terms appearing in the kernel (3.12), conservation entails three relations 
between them so that we need two independent functions. Furthermore, conservation should 
be automatically guaranteed by the form of the differential operators. We therefore consider 
the operators 

Sab = dadb-r]aba, (4.1) 

which fulfil d°'Sab = identically (i.e., acting on any tensor Tm—n)- An ansatz for the stress 
tensor two-point function which respects the symmetries is then given by 

{Tab{x)T,dm = '^{Sa(cSd)b - SabScd)f{x^) + SabScd9{x^) , (4.2) 

where now the functions / and g are unconstrained.'^ To determine /(a;^) and g{x^), we make 
for them the ansatz 

_ /(.-)n.-,--^ -'"-;>::'|:^,^'-' n.,^ ,4.3, 



^As a side note for later generalization to curved space backgrounds, we note that R — Sabh""^ up to linear 
order (see appendix A). 
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and analogue for g{x'^). Plugging this into the general ansatz (4.2) and comparing with the 
regularized result (3.11), we obtain five relations for F{z) and G{z) which are solved by 



F(z) 



G{z) 



2z — n 

2{z + 2- n){z + 1 - n){2z + 2 - n) 

(n - 2zf 
2{z + 2- n){z + 1 - n) ' 



(4.4) 



so that we obtain 



{Tab{x)Tcd{^)) = {Sa(^cSd)b - SabScd)x 



X J {m'Yix') 



2^..^2^.+2-n^(n - 2 - z)r (f - 1 - z) T (f + 1 - z) T{-z) dz 



222+47r"r(n + 2-2z) 27ri 

2v/.2,.+2-nr(n - 2 - z)r2 (§ + 1 - z) r{-z) dz 



22^+37r"r(n + 2-2z) 



2tt\ 



(4.5) 



The contour runs left of 'Rz = 0, where the first pole is encountered. In four dimensions, the 



C 



Figure 2. The first McUin-Barncs integral appearing in the stress tensor correlation function. There 
are simple poles dX z — k (black dots), atz = n— 2 + /c (white dots) and atz = ^ — 1 + (black 
squares) and double poles atz = ^ + l + fe (white squares), shown here for n = 3.75. The original 
contour C runs left of all poles. We lift it over the pole at z = to obtain the contour C* and the 
residue of this pole (shown in grey). In the remaining integral over C* , we can take the limit n 4, 
and the poles flow together (shaded circles). 



residue of this pole is proportional to (x^)^^ which is still not a well-defined distribution. 
However, the residues of all other poles are (with the proper Feynman or Wightman prescrip- 
tion) well-defined distributions in four dimensions. We therefore lift the contour over this 
pole (see figure 2) to extract the problematic term and take the limit n — )• 4 in the remaining 
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integral. This leaves us with 

(4.6) 

To treat the problematic term, we apply the reduction formula (2.5) to obtain 

i^'?-'' = ^( ^ U{xy-- , (4.7) 

l[n — 6)[n — 4) 

and the limit n ^ 4 could be taken to obtain a well-defined distribution, were it not for the 
explicit power (n — 4)~^. It is now that the difference between time-ordered functions and 
Wightman functions becomes important. Namely, for the (massless) Feynman propagator 
we have (equation (2.4) for a = 1) 

n{0\T4>{x)m\0) = □ . \ ^ = i5"(x) , (4.8) 

while for the Wightman function it results 

□ (O|0(x)<A(O)|O) = □ 7 ^ \ , -nZT = . (4.9) 

\ 47r2 (x"^ + lO sgn t) 2 ^1 

We therefore add to (4.7) zero in form of /x""^[(a;^)^~ 2 — (x^)^~2] with an arbitrary mass 
scale and obtain 



x2)2- = □ (^^)^-" _ ^i^-\x') 



„2\3-n ,,n-4/^2\l-5 



,71—4 



2(n-3)(n-4) ^ ^ v ; j 2(n-3)(n-4) 

1 /ln(//V)\ 2M " ... .^ 

□ ^^V^ + 7T7 ^^7^7 ^ □(x2)i-2 + 0(n - 4) , 



4 V 7 2(n-3)(n-4) 

which for the Feynman propagator reads 

„2,.nA2-n_ 1 ^ / M/X^^E^ + iO) \ , ^ , ^ 2 



(x^+iO)^- = -- □ J +ivr- ^_ _ 2 + 7 + In^ + 21n/.J <5"(x)+0(n - 4) , 

(4.11) 

while for the Wightman function we just have 

/ 9 s9 r, 1 ^ / ln(u^x^ -I- iO sgn f) \ , . 

{x^ + iO sgnt)^-" = -- U[ — — ^ + 0{n - 4) . (4.12) 



4 \ -|- iO sgnt 

Those formulas just effect the extension of the problematic term to the diagonal x^ = 0. 
That is, they coincide with the original term for all x'^ ^ 0, but are well-defined distributions 
(in four dimensions) also when smearing with test functions which have support at = 0. 
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The procedure described above is just dimensional regularization and renormalization in 
coordinate space [21, 31, 32], and the connection with the usual momentum space techniques 
is described in appendix B. 

Summarizing, we can decompose the stress tensor two-point function into a singular 
and regular part, where the singular part is local and only appears for the Feynman (and 
Dyson) propagator. The singular terms have to be subtracted with counterterms, while the 
regular part is the renormalized stress tensor correlation function. Schematically, we write 
for the various two-point functions (T~^ stands for anti-time ordering) 

{rTab{x)TUO)) = i{Tat,T,,r^ix) + {TM^'^ix' + iO) 
{Tab(x)T,dm = {TabTcdY'^ix^ + iO sgnt) 
{Tcd{0)Tab{x)) = {TabTcdVHx^ - iO sgnt) 
{T-'Tabix)TUO)) = -iiTabT^r^x) + {TabTcdr^ix^ - iO) 
where the regular part is given by 

{TabTcdV^^ix ) = {Sa(cSd)b - SabScd) ^ 



(4.13) 



□ 



V 38407r4x2 



+ 



+ SabScd 



2^., 2^.-2^(-^)^(l-z)^(2-z) dz 



{m'Yixy 



27ri 



/ln(/i^x^ 



/ / 2x2/ 2Nz-2r(--z)r(2 - 2;)r(3 - z) dz 
2 1+1 ("^ ) (x^)^ 2 \ ; \ 
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c* 



2567r2r (^ - z) 



and the singular part reads 



(TabTcd) 



sing/ 



1 



46 



9607r2 Vn-4 15 



+ 



1 



2407r2 



n 



47 
30 



7 - In TT 2 In /xj {Sa(cSd)b - SabScd) (^"(a:) 
-7-ln7r + 21n^^ SabScd S^'ix) . 



27ri 
(4.14) 



(4.15) 



We now show that those singular terms can be subtracted by known counterterms which are 
quadratic in curvature tensors. For this we need that up to surface terms and to linear order 
in the perturbation hab we have 



h''''SadSbch'''d''x 

h'^'SabScdh^'d^X 



R'''"'RabcdV^d''x = 4 / R'^'RabV^d'^X 



jab ] 



i?2 



-gd X 



R' 



-gd X, 



(4.16) 



which can be easily proven using the expansions given in appendix A. We therefore obtain 



4807r2 



h''\x){TabTcdr''^{x - x')h'''{x') d^xd^x 



6 I R^''RabV^gd''x + 



n 



15 



7 — In TT -|- 2 In /X 



Zrt itab \ -ft 



-g d^x . 
(4.17) 



The infinite part of this expression coincides with well-known results [33, 34]. 
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5 General Hadamard states 



A general state of a free quantum field is called Hadamard if near coincidence its two-point 
function is of the form 



G{x, x') 



u(x, x') 
47r2(x — x'Y 



+ x) ln((x — x) ) + w{x, x) , 



(5.1) 



where u{x, x'), v{x, x') and w{x, x') are smooth functions with u{x, x) = 1 [35, 36]. Hadamard 
states have a lot of nice physical properties (for example, the Hadamard form is preserved 
under Cauchy evolution, and vacuum states in static spacetimes are Hadamard), and it has 
been shown in examples that non-Hadamard states have unpleasant behaviour (such as an 
infinite stress tensor expectation value with usual renormalization procedures) [37-46]. In 
particular, the standard Minkowski vacuum state |0) is a Hadamard state, which follows by 
expanding the two-point function (2.7) in four dimensions near = 0, 



Go(x2) 



(27r)Vx2 



Kl I mv x^ I ~ 



1 m 

■ + ■ 



Ati'^x'^ 167r2 



l-27-ln(m^x74))+0(x2) . (5.2) 



Since v{x) is fully determined by the geometry of the underlying spacetime, two Hadamard 
states differ only by a smooth function. Especially, the two-point function in any Hadamard 
state |H) can be written as 



Gh(x, x') = Go((x - x'Y) + WYi{x, x') 



(5.3) 



where W}i{x,x') is a smooth function. If we therefore calculate the stress tensor correlation 
function in a general Hadamard state, we obtain (using (3.7) and (5Go)wh = diG^wii) — 
Goidwu)) 



{Tab{x)Tcd{x')) ^ = {Tab{x)Tcd{x')) ^ + -m'^T]abr]cdWHWu + m'^rjab'ncdGoWu 

- ^^Wll [(d^pwn) {d^qwn) + 25f (Cod^^ wn) - 2Go«w;h 

- VcdmS^l [(5>h) {dfwn) + 2dt {G,dlwn) - 2G,dtdlwu\ 

+ ^Vabra [dtd^ {Godfd^^'wn)+Godtdfd;'d^^wn 

- ^VairfJ (Godfd^'cf^wn) + d^' {Godldfd^^w^ 



with the short-hand notation 



Vab 



1 



-.VabV 



kl 



(5.4) 



(5.5) 



and wh and Go depend on x and x' . If we replace {Tab{x)Tcd{x')) q by its renormalized 
value (4.14), this is already a well-defined distribution, because derivatives of a smooth 
function (of wu) are still smooth and a distribution (Go) can be multiplied by a smooth 
function to yield another well-defined distribution. Especially, the counterterms (4.17) which 
one needs to subtract the divergences as n — t- 4 are the same for all Hadamard states. 
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As a simple example, with the scalar field decomposed in the standard creation and 
annihilation operators [30, 47], we consider the states \a, k) = exp[ia(a^ + afe)]|0) for real a 
and in the massless case m = 0. Using the Baker-Campbell-HausdorfF formula [48], we see 
that they are normalized, and an exercise in commutation relations^ then gives the two-point 
function in such a state which is 

ia^k\(j){x)(t)ix)\a,k) = (O|(/)(x)0(x')|O) + — j- [sin(A;x) sin(A;3;')]fco=|fc| (5-6) 

|fe| ' ' 

The smooth difference between the Minkowski vacuum two-point function and this special 
Hadamard state is therefore given by 

Wa,k{x, x) = —— sva.{kx — \k\t) sin(fca;' — \k\t') , (5-7) 
and its renormalized stress tensor two-point function reads 

{TMTM))Z = - iSa,.S,,, + 35.,5..) □ [^0^^, 

+ 8a'^^^^^cos\kx)cos\kx') (5.8) 
+ -^TTT cos(A;x) cos{kx')r]'^l'i]l^kkkp didq{x - x')~^ , 

where k^ = \k\\s understood. 

6 Comparison with earlier results 

In earlier calculations, only partial results were obtained. The general regularized case (in- 
cluding a general coupling to the curvature ^Rcf)^ in the action) was considered in [49], which 
coincides with our result (3.9) if one replaces G(x^) with the explicit form (2.7). However, 
no renormalization was performed. 

In [34, 50] renormalization was performed, but only for the massless case. In this limit, 
the Mellin-Barnes integral in the renormalized result (4.14) vanishes because the modified 
contour C* has > 0, so that we are left with the first term only 

{TabTcdy'^ix^) ^ -38^(^a(c5d)6 + □ {^^^^^ (m ^ 0) . (6.1) 

By inserting the explicit definition of the operators Sab (4-1), one obtains exactly the re- 
sult in [34], whereas to show agreement with the result of [50] we must extract another 
d'Alembertian operator 

. „ . 8 [In(/i2x2) + 1] 
□ \n^U^x^) = ^ J ^ . 6.2 

Since there only the Wightman function was considered where no additional local terms arise, 
we obtain 

□ □ In^ {^?x^ + iO sgn t) = 8 □ f (g.g) 

\ + iO sgnt / 



^Derive first [{alj + aq)^ ,<f){x)] = — ifc(2/|q|)2 [sin(ga::)]qO^|q[(aq + ciq)'° ^ and then commute the exponential 
with the fields. 
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and find full agreement with our result. 

The only fully renormalized result for general mass and coupling to curvature was de- 
rived in [12], but in momentum space where the calculation is cumbersome as stated by the 
authors themselves. To compare this with our coordinate space result, we need to take the 
Fourier transform for all the different prescriptions. Since the integration over the contour C* 
in the finite part (4.14) is absolutely convergent, we can exchange this integration with the 
Fourier transformation and Fourier transform the integrand. For the Feynman prescription 
we may use the formula (2.4), while for the Wightman prescription formula (2.10) is the right 
one. To fully transform the renormalized result (4.14), we need also the Fourier transforms 
of ln(x^)/x^ for the different prescriptions. This can be obtained by writing 



ln(x^ 



lim - 

e-5-0 e 



1 



1 



2\l+e 



(X2) 



2\l+2e 



and by using (2.4) and (2.10) we calculate 
ln(x2 + iO) 



{x^ + iO) 
ln(j;^ + iO sgnt) 



4i7r^ 



(p2 - iO) 



(X2) 



[27-ln4 + ln(p2 -iO)] 



(6.4) 



g-ipx d^^ = 87r3G(p°)G( V 



(6.5) 



(x^ + iO sgn t) ~ ' ' ~ ' 

Combining all of the above, for the Feynman prescription we obtain from (4.14) 



/ 



cd. 



reg, 



+ iO)e-'P^ d^X 



+ 



1 

5127ri 
i 

2567rt 

where the operators Pat = 
and we defined 



^^(P-(cPd)b + 3PabPcd) (27 - ln4 + ln(/ - iO) - ln{fi^)) 

r(z)r(-z)r(i - z) dz 



{Pa{cPd)b 



PabPcd) / (T 



(6.6) 



-PabPcd 



a' 



T(z)r(-z)r(3- 
ra-z) 



z) dz 
27ri 



27ri 



PaPb are the Fourier transforms of the operators Sab (4.1) 
4m2 



0- = . (6.7) 

p2 - iO ^ ^ 

Note that because of the product of F functions T(z)T{—z) in the numerator, the standard 
Mellin-Barnes contour C is not defined for those integrals, but the modified contour C* (see 
figure 2) works. We may then close the contour to the left if |cj| > 1 or to the right if lo"! < 1 
and sum the residues to obtain 



(Ta6red)'"^^(x^ + iO)e-'f^d"x 



+ 



+ 



4807r2 



2407r2 



(PaicPc 



{c^d)b 



PnhPr 



ab^cd) 



46 
30 



+ 7 + ln 



m 
2^ 



+ 



1^F(|)F(A; 



7^ 

fc=0 



PabPcd 



47 
30 



f m\ 1 \ - 

+ 27 + 21n(-j 



r(i) 



T[k + l) 



-(-1)V 



k-k-l 



(6. 
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With hindsight, we rewrite this as a sum over 1 + a using 

fc=0 m=0 ■ \k=0 ' ) 

The sums in the expression (6.8) then can be done in terms of GauB' hypergeometric function 
which reduces to even simpler functions for the cases considered, and it results 



v r(i)r(A: + i) (i + oT r(|)r(A. + m + i) 

fc=0 1^^+ 2J „^=o Vfc=0 " V'^ ' 1> 

7 

~ ~15' 

(6.10) 

and 



[3 + 5(1 + a) + 15(1 + o-)2] + 7(1 + (t)^ artanh \/l + a 



y r(i) r(fc + 4) fc^-fc-i ^ _1(94 _ + 45^2^ 

+ - \/l + cj(8 - 4cr + 3cr^) artanh Vl + • 

For the Fourier transform of the Feynman prescription stress tensor correlation function (6.8) 
we therefore obtain 

+ -^:^{Pa{cPd)h - PabPcd) 7 + In - ^ct{7 + 3a) + (1 + a)i artanh Vl + ^ 



+ 9607r2^'^'^'='' 



87 + 8 In + 3ct(1 - cr) + Vl + cr(8 - 4cr + 3cr2) artanh ^1 + o" 



(6.12) 



The result of [12] is given for the minimally coupled case and for the Feynman prescrip- 
tion as (expressed in our notation) 



(r„,r,,rs(x2 + io) = ^ 



^ (35,(,5,), - SatS^d) + ^)'0(/)e'^^^ 



d^p 



4 



+ ^SabScd I {a- 2)20(/)e 

"T-^ ^TTTT (35a(c5'rf)fc - SabScd) + -^SabScd^ J e^^^ — 



135 ' ^^^ao^caj J ^2 (2^)4 

(6.13) 

up to local terms, which anyway depend on the exact renormalization scheme used. In their 
appendix we find 

</,(/) = -2 + TTTa In f^^^S^^ = -2 + 2Vl + cr artanh Vl + a, (6.14) 

VVl + cr - ly 
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and so their result for the Feynman prescription is (up to local terms) 

+ Y^iPa{cPd)b - PabPcd) -^(J{7 + 3a) + (1 + (j)t artanh ^/^^ 
+ ^^^P^bPcd [3ct(1 - cj) + (8 - 4cr + 3cr2)Vl + cT artanh Vl + o- 



(6.15) 



Up to a global factor of 4, this coincides exactly with our result (6.12). 
For the Wightman prescription, we do the same steps to obtain 



,5 



I {TabTcdY'^'ix^ + iO sgnt)e-'f- d^x = ^(P,(,P,), - p^,p^,)e(/)e( V)(l 

(6.16) 



+ 7^^abPcde(p°)e(V)\/r+^(8 - 4cT + 3^2) 



For this prescription, the result of [12] reads (expressed in our notation) 



(T„6Te,)'-^S(^2 ^ -0 gg^^) ^ ^ (35„(,5,), - SabScd) + a)^I{p)^^-^, 

+ "^SabScd I {a- 2)2/(p)e'P-'^ 



(6.17) 



(2^)4 ' 



where 



i{p) = ^(^(1 - sgnpO)e( V - 4m2)7rT^ = -i^e(-/)G( V)e(i + a)Vl + ^, 

(6.18) 

and therefore they have 

J {TabTcd^'^ix^ + iO sgnt)e-'P^ d'x = 

+ ^(PaicPd)b - Pa6^'cd)e(-pO)e( V)e(l + a)il + a)l (6.19) 
+ 7^PabPcd@{-p'')@{-p^)&il + a){8-Aa + 3a^)Vl + ^ , 

in full accordance with our result (6.16) except for the sign of p*^ (which probably is only a 
typographic error) and the additional Heaviside function 0(1 + a). 

In summary, we may conclude that up to minor differences our result is in full agreement 
with earlier findings. 



7 Conclusions 



We have presented an explicit procedure to obtain renormalized two-point functions of stress 
tensors directly in coordinate space, and have given the full renormalized result (4.13) for 
the case of a minimally coupled scalar field in Minkowski spacetime. Our result provides 
an important check on earlier partial results and results obtained in momentum space. This 
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result can now be used to determine quantum matter corrections to the background geometry 
which give rise to a variety of interesting effects, also in curved spacetime [51-71]. 

Our method rests on two points: first, a linearization formula (2.12) for the product 
of two propagators and their derivatives, representing them as a Mellin-Barnes integral; and 
second, the extraction of differential operators to bring the correlation function into a form 
well-suited for renormalization (4.5). The method works directly in coordinate space and 
so should be more easily generalizable to curved spacetimes where no Fourier transform is 
defined, which is especially important in the cosmological context. In fact, for the relevant 
case of a de Sitter background, a similar linearization formula for the product of two (undif- 
ferentiated) propagators has been derived [72-74] and used in some applications [75, 76] (see 
also [77, 78]). Of course, the linearization formula derived in this article depends on the fact 
that the two-point function of a scalar field state invariant under the Poincare group only 
depends on the invariant distance (x — x')^. For a generalization to curved spacetimes, this 
will not the case; however, say for a FLRW universe with scale factor a{r]) depending on the 
conformal time r/, is is conceivable that a similar Mellin-Barnes representation exists with the 
kernel additionally depending on r] and rj' . For such a representation, presumably a lineariza- 
tion formula can be found and the short-distance singularity extracted as in Minkowski space 
(since the kernel, even though it depends explicitly on time, has a finite limit as x' — ?■ x). 

The second important point - extracting differential operators - exploits the fact that 
conservation restricts the form that the stress tensor two-point function can take. For a 
Poincare invariant state, taking into account conservation reduces the number of independent 
functions from five to two, and we extract those operators such that conservation is fulfilled 
automatically. The remaining two functions are then in principle completely arbitrary, so 
that one can perform renormalization separately on the most divergent terms which are the 
only ones that need to be renormalized. It has the further advantage that the strength 
of the divergence in those functions is reduced to the one of a product of undifferentiated 
propagators, which need counterterms that diverge like (n — 4)~^ as appropriate for a one- 
loop calculation. For the stress tensor correlation function of the massless, minimally coupled 
scalar, this idea is used in [34, 50] in Minkowski space and in [79] in de Sitter space; it also 
arises quite naturally in Fourier space in a Minkowski background [12]. However, in curved 
spacetime it is not easy to find the appropriate generalization for the operators Sab (4-1), 
whose definition was motivated by their symmetry and vanishing divergence. As noted 
before. Sab with both indices referring to the same point can presumably be substituted 
by the expansion of the Ricci scalar up to linear order in the metric perturbations R = 
Sabh!^^ . However, the term Sa(cSd)b where the index pairs {ah) and {cd) refer to different 
spacetime points will generally fail to be divergence-free when acting on a scalar function 
when generalized in the obvious way by substituting parallel propagators for the metric; 
it is only in Minkowski space where the parallel propagator coincides with the metric that 
its divergence vanishes. In [79], the differential operator obtained by the expansion of the 
Weyl tensor was proposed but yields unwieldy results.^ From the gauge invariance of the 
interaction habT"'^ and the fact that the curvature tensors (minus their background value) 
are gauge-invariant, one may strongly suspect that those differential operators have to come 
from the expansion of curvature tensors, so that maybe the choice of the Riemann tensor 
instead of the Weyl tensor will give nicer results. 

^This includes strong divergences like (n — 4)""^ in intermediate steps, which however may be an artefact 
of the concrete calculation — the final result is the expected one, and the counterterms coincide with well 
established ones. 
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The calculation of the stress tensor correlation function for other types of fields can 
be done with the same ease; since their propagators are obtained by taking derivatives of 
the scalar propagator, the exactly same Mellin-Barnes representation can be used. The only 
nontrivial issue (in the author's opinion) is to check if the integrand in the analogue of 
equation (3.11) is conserved on its own so that the rest of the procedure goes through. 
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A Curvature tensor expansions 

In this appendix we present the expansion of curvature tensors around a Minkowski back- 
ground up to linear order in the perturbation. 

The perturbed metric and the Christoffel symbols are given by 

gab = IJab + hab 
gab ^ ^ab _ f^ab 

1. (A.l) 



v-g = 1 + -h 

where the indices on hat are raised using the background Minkowski metric, and the curvature 
tensors follow straightforwardly 

Rab = \ {dad%d + dbd^'had " daObh - Dhab) (^-2) 

R = dadbh"^ - ah . 



B Momentum space regularization and renormalization 

In this appendix we show how dimensional regularization and renormalization in position 
space is related to the usual momentum space techniques. For simplicity, we will treat the 
massless case, which also is the only term we need to renormalize after extracting differential 
operators from the stress tensor two-point function. Using the formula (2.4) for the Fourier 
transform of the massless Feynman propagator, we obtain the well-known expression 

1 47rt 1 

^e-'P^d"x = -i— . (B.l) 

(x2 + i0)t-i r(5-l)p2_iO ^ ^ 
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To the product in coordinate space corresponds a convolution in momentum space 

(x2+i0)t-l (x2 + i0)t-l'' ^ = -r2(|-l)y (p-<;)2_i0g2_i0^ 9 

. 4^-t7rtr(2-§) 1 ^^-^^ 

' r(n - 2) (p2 _ io)2-t ' 

which can be done in n dimensions but is divergent as n — )• 4. In momentum space one now 
expands around n = 4 to get 

/ (^2 + -o)n-2 ^"^^^ d"x = ivr^ - 2 + 37 - ln4 + In^r + ln(p2 - iO)) + 0{n - 4) . 

(B.3) 

The Fourier transform of the logarithm can be taken using the formula (2.4) by rewriting it 
as 



(B.4) 



_ /27-ln4 + ln(x2 +iO)\ 
~ V 47r2(x2 + i0) ) ' 

and by using □[47r2(2;2 + iO)]^^ = 15^ (x), we get combining all of the above 

- '-^ - 2 + ^ + 1"') ^*(-) - i ° + o(" - ^) ■ (B-^) 

This is of course exactly the result we obtained previously in coordinate space (4.11). 
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